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Abstract The CMB angular temperature fluctuations observed by Cobe and Wmap enable us to place a lower 
limit on the spectral distortions of the CMB at any angular scale. These distortions are connected with the simple 
fact, that the superposition of blackbodies with different temperatures in general is not a blackbody. We show, 
that in the limit of small temperature fluctuations, the superposition of blackbodies leads to a y-type spectral 
distortion. It is known, that the CMB dipole induces a y-type spectral distortion with quadrupole and monopole 
angular distribution leading to a corresponding whole sky y-parameter of yd ~ 2.6 • 10~^. We show here that 
taking the difference of the CMB signal in the direction of the maximum and minimum of the CMB dipole due to 
the superposition of two blackbodies leads to a spectral distortion with i/opt — 12 yd — S.l ■ 10~^. The amplitude 
of this distortion can be calculated to the same precision as the CMB dipole, i.e. 0.3% today. Therefore it may 
be used as a source with brightness of several or tens of to cross calibrate and cahbrate different frequency 
channels of CMB surveys with precision of a few tens or hundreds of nK. Furthermore, we show in this work 
that primordial anisotropies for multipoles 2 < I < 1000 also lead to spectral distortions but with much smaller 
y-parameter, i.e. y ~ 10~^^ — 10~^. 
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1. Introduction 

The Wmap spacecraft measured the amplitude of the an- 
gular fluctuations of the cosmic microwave background 
(CMB) temperature with extremely high precision on a 
very broad range of angular scales, from 12' up to the 
whole sky. These temperature anisotropies and the exis- 
tence o£_th^j£Oustic_2eaks_wer^^ 
ago l|Peebles fc YuL Il970l: ISunvaev fc ZeldovichL Il970a|) . 
but only now after Boomerang, Maxima, Archeops, 
Wmap and many ground based experiments like Cbi, 
AcBAR, VSA, etc. we know their precise characteristic an- 
gular scales and amplitudes. With no doubt, only tiny 
fluctuations of the radiation temperature field have been 
observed, but nowadays with a precision to better than 
1% down to degree angular scales. 

It is commonly assumed that the spectrum in one di- 
rection of the sky is planckian and that only the tempera- 
ture changes from point to point. This follows from the na- 
ture of the main effects leading to the app earance of these 
fluctuations, i.e. the Sachs- Wolfe-effect llSachs fc Wolfd . 
Il967fl and the Doppler effect due to Thomson scatter- 
ing off moving electrons ( Sunvaev & Zcldovich, 197Da.) at 
redshift ~ 1100. However, as will be demonstrated below, 
there are spectral distortions in second order of AT/T. 
These distortions are inevitable when the CMB is ob- 
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served with finite angular resolution or when regions on 
the sky containing blackbodies with different tempera- 
tures are compared with each other. 

The CMB missions mentioned above have shown that 
there are fluctuations of the radiation temperature on the 
level of AT ~ /iK— mK over a broad range of angular 
scales. One may distinguish two basic observational strate- 
gies: (i) absolute measurements, where the beam flux in 
some direction on the sky is compared to an internal cali- 
brator (Cobe/Firas) and (ii) differential measurements, 
where the beam flux in one direction on the sky is com- 
pared to the beam flux in another direction (Cobe/Dmr 
or Wmap). In the flrst strategy one observes a sum of 
blackbodies (SB) due to the average over the beam tem- 
perature distribution, whereas in the second two sums of 
blackbodies are compared with each other. Under these 
circumstances we will in general speak about the superpo- 
sition of blackbodies, i.e. the sum and difference of black- 
bodies with different temperatures. 

Any experiment trying to extent the great success of 
the C obe/Firas instrum e nt, which placed strict upper 
limits llFixsen et all Il996t IPixseii fc Matheil l2002|) on a 

possible n- ("Sunvaev fc ZeldovichMl970lJi), < 9 • 10-^ 
and y-typc (Zcldovich & Sunvaev, 1969), |y| < 1.5 • 10-^ 
CMB spectral distortion, will only have a flnite angular 
resolution and would therefore observe a superposition of 
several Planck spectra with different temperatures corre- 



spondmg to the maxima and mmima on the CMB sky as 

measured with Wmap. 

It is known l|Zeldovich. Illarionov fc Sunvaevi Il972l) . 

that in the case of a gaussian temperature distribution this 
win lead to a spectral distortion indistinguishable from 
a y-distortion, with corresponding y-parameter which is 
proportional to the dispersion of the temperature distri- 
bution. Since the temperature fluctuations of the CMB 
indeed are gaussian, this implies, that the corresponding 
spectral distortions averaged over large pieces of the sky 
should be of y-type. But here we are interested in the case 
of measurements with angular resolution of a few arcmin- 
utes to degrees. In this situation, we deal with the limited 
statistics of finite regions with different mean tempera- 
tures and therefore it is not obvious, what type of spec- 
tral distortion would be induced in each small patch of 
the sky. As mentioned above, a similar situation arises 
when we compare the signals from two regions on the sky, 
i.e. take the difference of the intensities as is usually done 
in differential observations. Below it will be shown, that 
for any observation of the CMB temperature fluctuations, 
there will be unavoidable spectral distortions due to the 
difference in the temperature of the radiation we measure 
and compare and that these distortions will be indistin- 
guishable from a y-type- distortion. The biggest distortions 
arise due to the CMB dipole. 

The unprecedented high sensitivity of future or pro- 
posed space missions like Planck and Cmbpol or ground 
based instruments under construction like Apex, the 
South Pole Telescope (Spt), the Atacama Cosmology 
Telescope (Act) and Quest will offer ways to investi- 
gate tiny secondary CMB angular and spectral fluctua- 
tions and should therefore add a lot to the success of pre- 
vious missions. One target will be the measurement of 
the SZ effect from clusters, proto clusters or groups o f 
galaxies and superclusters l|Sunvaev fc ZeldovichI Il972^ 
or signatures from the first sta rs in the universe 
l)Oh. Coorav fc KamionkowskiL l20oi) . In the future CMB 
experiments will be so sensitive that it will be possible to 
investigate in detail the imprints of reionization and the 
tra ces of energy relea se in the early universe. 

iBasu et al.l l)2003l) proposed a method to constrain the 
ionization history of the universe and the history of heavy 
element production using the properties of resonant scat- 
tering of CMB photons in the fine structure lines of oxy- 
gen, carbon and nitrogen atoms and ions produced by the 
first generation of stars. The strong frequency dependence 
of this effect permits to extract the undisturbed angular 
dependence of the frequency independent primary tem- 
perature fluctuations and thereby avoid cosmic variance. 
By comparing the signals in different frequency channels 
it is possible to investigate the contributions of the lines 
of different ions at different redshifts and therefore to ex- 
amine different scenarios of element production and ion- 
ization histories in the low density regions of the universe, 
with overdensities less than <^ 10"*. The sensitivities of 
Planck and Act should be sufficient to detect the sig- 
nals imprinted by the effects of resonant scattering, but 
the crucial point for the successful measurement of any 
small frequency dependent signal is the cross calibration 
of the different frequency channels down to the limits 
set by the sensitivity of the experiments. Full sky mis- 



sions like (JOBE/Umr or Wmap normaUy use the CMB 
dipole and its annual modulation to check the calibra- 
tion of their instruments down to a level of /iK, whereas 
experiments with partial sky coverage like Boomerang 
are directly using the CM B dipole for calibration issues 
l|de Bernardis et all l200o|) . But both methods permit to 
cross calibrate different frequency channels only in first 
approximation assuming that the dipole has a planck- 
ian spectrum and the same amplitude at all frequencies. 
Unfortunately, this precision of the cross calibration will 
not be sufficien t to detect the sig nals from the dark ages 
as discussed bv iBasu et all l|2003|) . 

It is known, that the motion system relative to 
the CMB restframe in addition to the dipole gener- 
ates (in second order of v/c) a small monopole and 
quadrupole contribution to the CM B brightnes s of the sky 
in th e restframe of the o b server JSunvaev fc ZeldovichI 
198(i Ide Bernardis et ail llQQOt iBottani et all Il992|) . 



Sunvaev fc ZeldovichI l|l98fJ) . when they were discussing 



the radiation field inside a cluster of galaxies mov- 
ing relative to the CMB restframe, have shown that 
the corresponding dipole induced quadrupole has a 
non planckian spectrum , which was then later derived 
bv ISazonov fc SunvaevI l|l999f) . iKamionkowski fc Knoxl 
'20039 later applied this solution to the case of our mo- 
tion relative to the CMB restframe and proposed to use 
the dipole induced quadrupole for calibration purposes. 

The solution of ISazonov fc SunvaevI l|l999l ) is valid in 
the case of a narrow beam observations. In this paper we 
choose an independent approach, which is based on the 
superposition of blackbody spectra with different temper- 
atures, to look for the maximal and minimal spectral dis- 
tortion obtainable from CMB maps. Our method allows 
us to calculate the value of the y-parameter for the dipole 
induced monopole and quadrupole for a beam with finite 
width or equivalently for any average of the signal over ex- 
tended regions on the sky. Most importantly we show that 
the difference of the sky brightness in the direction of the 
maximum and minimum of the CMB dipole, correspond- 
ing to the maximal difference of the radiation temperature 
on the CMB sky, leads to a y-type spectral distortion with 
associated y-parameter of ?/opt = 3.1 • 10^^. We propose 
here to use this spectral distortion arising due to the CMB 
dipole to cross calibrate the frequency channels of a CMB 
experiment in principle down to the level of a few tens of 
nK. We discuss different observing strategies in order to 
maximize the inferred spectral distortion (Sect.0l- 

In this paper, we first give a short summary of the 
basic equations necessary in the following derivations and 
define some of the terminology used (Sect. O. We then 
discuss the underlying theory for small spectral distor- 
tions (Sect.O and show that in this limit even the distor- 
tions arising due to the superposition of two blackbodies 
(Sect. 01 with close temperatures are well described by a y- 
type solution. Furthermore, we discuss in detail the spec- 
tral distortion due to the superposition of Planck spec- 
tra with different temperatures arising from to the CMB 
dipole (Sect. jS]) and from the higher multipoles (Sect. O 
using generated CMB sky maps for the Wmap best fit 
model. We discuss the spectral distortions arising in dif- 
ferential measurements of the CMB temperature fluctua- 
tion (Sect.[7|) and ways to use the spectral distortions in- 



duced by the CM_B dipoie to cross calibrate the trequency 
channels of CMB experiments (Sect.lS)). We end this work 
with a discussion of the consequences of the obtained re- 
sults for some of the highly demanding tasks which may 
be addressed by future CMB projects (Sect.EJ and finally 
conclude in Sect. 1101 



2. Basic ingredients 

2.1. Compton y-distortion 

In the non relativistic limit, the comptonization of the 
CMB photons by hot, isotropic, thermal electrons with 
Compton y-parameter 



2 o'T rie 



dl. 



(1) 



where To is the temperature of the electron gas, ax is the 
Thomson cross section and n^. is the electron number den- 
sity, leads to a y-distortion ( Zcldovich fc Sunvacv.. .1969.) : 
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for y <C 1. Here x = hv/kTo is the dimensionless fre- 
quency, v is the photon frequency, Tq is the temperature 
of the incoming radiation and A/ — I — Iq denotes the dif- 
ference between the observed intensity / and undisturbed 
CMB blackbody spectrum Iq — B,y{To), with 



BAT) 
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(3) 



In the Rayleigh Jeans (RJ) limit, i.e. a; <C 1, equation 
(0) simplifies to AT/TqIrj = —2y. Since in this limit 
a blackbody of temperature T is approximately given 
by B^Yi} = ^v'^kT^ it is convenient to compare the 
comptonized spectrum to a blackbody of temperature 
Trj = To [1^2 y], making the relative difference (0 vanish 
at small v. 
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Here we now defined x — hv/kT^:^ and / = B,j{T^j) is a 
blackbody of temperature Trj. The corresponding differ- 
ence in the radiation temperature can be written as 



AT 



y-gy{x) = y 
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e* - 1 



(5) 



This equation shows, that the most important character- 
istic of a y-distortion is given by the function 
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(6) 



Whenever the function g{x) characterizes the frequency 
dependent part of the relative temperature difference we 
may speak about a y-type spectral distortion. 



z.Z nelation between temperature and intensity 

To obtain the relative difference in temperature © corre- 
sponding to a y-distortion from the relative difference in 
intensity I0J the relation 
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AI dlnT 
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was implicitly used. In CMB measurement this relation 
is commonly applied to relate intensity differences to tem- 
perature differences, assuming that the difference in inten- 
sity obeys ^ 1 or equivalently <^ 1. The sensitiv- 
ity of future CMB experiments will be extremely high in 
a broad range of frequencies and angular scales. Therefore 
it is important to understand what corrections arise in 
next order of and and what kind of distortions are 
introduced by using relation . For this we consider the 
simplest case, when we are comparing two pure blackbod- 
ies with different temperatures T and T' and correspond- 
ing intensities / = B,,{T) and /' = B^{T'). Then using 
(PI the relative difference in intensities can be written as 
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where we defined x = hv/kT and A = (T' - T)/T. 

Now, using Taylor expansion up to second order in A 
from equation ||HJ) we obtain 



A/ 



X e 



+ 0(A3), (9) 



where gy (x) is defined by equation Q . Now we apply the 
relation lO to infer the relative temperature difference 



AT(f) 
T 



= A- 



9yix) ■ — 



0{A^). 



(10) 



This result clearly shows that in second order of A a y- 
distortion with y-parameter y — A^/2 is introduced. This 
distortion vanishes at low frequencies (x —^ 0) and in- 
creases like gy{x) ^ x in the Wien region (x ^ 1). Writing 
A + gy{x) ■ A'^/2 = A-[l+gy{x)- A/2] this implies that for 
any given A there is a critical frequency Xy which fulfills 



9y\Xy) — 



(11) 



and above which the contribution of the y-distortion rel- 
ative to the first order term becomes bigger than e. The 
value of e is directly related to the sensitivity of the exper- 
iments. In Fig.njaij, has been calculated for different e. At 
frequencies lower than Xy the relation IJJ) is exact within 
the sensitivity of the experiment, whereas for x > Xy the 
second order correction has to be taken into account. In 
the regime ^ 4 for estimates the simple approximation 
Xy =2 e/|A| -I- 2 may be used. 

To give an example, the biggest temperature fluctua- 
tion on the whole CMB sky is due to the CMB dipoie. 
Comparing the maximum and the minimum of the dipoie 
corresponds to A ^ 2.48 • 10""^. If some CMB experiment 
is able to accurately measure temperature fluctuations on 
the level of fiK, contributions of e ~ 10"''^ relative to the 
dipoie signal can be distinguished leading to Xy^d ^ 1- This 
estimate shows, that at high frequencies spectral distor- 
tions introduced by the usage of the relation JT)) for the 




j.l. bum ot blackbodies 



Figure 1. function of A for different values of 

e according to equation ((TT|l : For e = 10~^ (solid), e = 
10"'' (dashed), e = IQ-^ (dotted) and e = IQ-^ (dashed- 
dotted). The right ordinate corresponds to Vy — kTo ■ 
Xy/h, with To = 2.725 K. 



dipole anisotropy should be taken into account in future 
CMB missions like Planck and Act. In Sect. we will 
discuss these distortions due to the dipole in more detail. 

Let us note here, that the temperature difference 
should not be larger than a few percent of the tempera- 
ture of the reference blackbody. Otherwise corrections due 
to higher orders in A will become important and lead to 
additional distortions (see Sect.|2|). Fortunately, the tem- 
perature differences on the CMB sky are sufficiently small 
to neglect these corrections. 

Equation H10(l also shows, that in general the inferred 
temperature difference is frequency dependent. At fre- 
quencies below Xy the inferred temperature difference is 
close to the true temperature difference A and frequency 
independent within the sensitivity of the experiment. In 
this case, we will speak about a temperature distortion or 
fluctuation, emphasizing that it is frequency independent. 
It is possible to eliminate the temperature distortion us- 
ing multifrequency measurements. For x > Xy frequency 
dependent terms become important, which we will hence- 
forth call spectral distortions. 



3. Small spectral distortions due to the 
superposition of blackbodies 

When the CMB sky is observed with finite angular res- 
olution or equivalently if the brightness of parts of the 
sky (not necessarily connected) are averaged one is deal- 
ing with the sum and more generally with the superposi- 
tion of blackbodies. Here we develop a general formalism 
to calculate the spectral distortions arising for arbitrary 
temperature distribution functions in the limit of small 
temperature fluctuations and derive criteria for the appli- 
cability of this approximation. We first discuss the basic 
equations necessary to describe the spectrum of the sum of 
blackbodies (SB) as compared to some arbitrary reference 
blackbody (Sect.|OJ and then generalize these results to 
the superposition of blackbodies (Sect.|^J). 



Follow ing the paper of IZeldovich. lUarionov & SunvaevI 
l|l972l) we express the total spectrum resulting for the SB 
with different temperatures T as: 

/(z.) = J R{T)B,{T)dT. (12) 

Here R{T) denotes the normahzed (/i?dT = 1, i? > 0) 
temperature distribution function, which will be used be- 
low to model the beam of some CMB experiment, and 
(T) is a Planck spectrum of temperature T as given by 
equation (PJ- Now we want to compare the spectrum /(j/) 
to a reference blackbody of temperature T^cf- Defining 
5 — {T — Ti-c[)/Ti-e[ and inserting T = T^d (1 + ^) into 
equation may be rewritten as 
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eS/(i-H<5) _ 1 



(13) 



where x = hv/kT^-ct is the dimensionless frequency and 
A — Sirk (? . CMB temperature fluctuations are of the 
order /^T/T ^ 10^^ — 10"'^. Therefore we are interested in 
the case when \5\ <C 1. This allows us to perform a Taylor 
expansion of equation H13|l : 
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5=0 



where here ^ = x/{\ + 5) was defined. Keeping only terms 
up to fourth order in 6 this simplifies to: 



B^{5)=I,,i ^gfc(^)5^ 



(15) 



fc=0 



Here we used the abbreviation /rcf = B^{5 = 0) for the 
reference blackbody spectrum of temperature Tj-ef and de- 
fined the functions gk{i) as 



go{x) = 1 
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(16a) 
(16b) 

(16c) 
(16d) 
(16e) 



where gy (x) is defined by equation Q • In Fig. [21 the fre- 
quency dependence of the functions gi{x) is illustrated. 
Deviations from a blackbody spectrum become most im- 
portant in the Wien region and vanish in the RJ region of 
the CMB spectrum. 
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Figure 2. Frequency dependence of gi/gi = /STi/T-R j 5^: 
Dotted line ATi/TrjJ, which is equivalent to a temper- 
ature distortion, solid line gy{x) — 2AT2/Tf{j(5^, which 
is equivalent to a y-distortion AT/yTo as given by equa- 
tion lO, dashed line AT^/Tjij and dashed-dotted line 
AT4/TRJ (5^. The function gy{x) becomes larger than unity 
for V > 146.2 GHz or a; > 2.57, larger than 5 for 
V > 396.7 GHz or i > 6.99 and larger than 10 for 
ly > 681.3 GHz or x > 12. For v = 100 GHz or x = 1.76 
it follows gy = 0.5. Also shown as vertical lines are the 
Planck LFI (dashed) and HFI (dash-dotted) frequency 
channels. 

Inserting equation H15|l into 1)12(1 the relative difference 
Al/Jrcf = {I — Ii-cf) / Irci between the SB spectrum and the 
reference blackbody /j-cf can be derived 
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9k{x) {S") , 



(17) 



where we introduced the abbreviation ((5*^) = J RS'^dT 
for the fc-th moment of the temperature distribution. To 
find the corresponding difference in temperatures from 
(|17l) one may use the relation As mentioned before 
the inferred temperature difference in general will be fre- 
quency dependent. 

In equation ((17|l the term proportional to the first mo- 
ment, {S^), corresponds to a temperature distortion which 
introduces a frequency independent shift in the measured 
temperature difference. It is possible to eliminate this con- 
tribution by multifrequency measurements, since it does 
not change with frequency. The term proportional to the 
second moment is indistinguishable from a Compton y- 
distortion as given by equation with y-parameter 



ys 



(18) 



Higher moments lead to additional spectral distortions 
making the total spectral distortion differ from a pure y- 
distortion. As will be shown below for the CMB these 
higher order corrections can be neglected. 



Limiting cases tor 

In the RJ limit {x ^ 1), starting from (|12(l and using 
relation ((Jj) one can find 
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I 



12 



(19a) 
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This result shows, that for i; — > the functions 52, .93, and 



94 all vanish like 



reflecting the fact that the sum of 



RJ spectra is again a RJ spectrum and that no spectral 
distortions are expected for 2; ^ 1. In the first step we 
have only assumed that a; ^ 1. Therefore, up to second 
order in x equation H19a|) describes the SB for any 6 to 
very high accuracy. For S ^ I equation (|19a|l takes the 
form AT/TicflRj ~ (S) ■ [1 + j^]. For the approximation 
(|19b(l we have in addition to a; <C 1 assumed that 5^1. 
In this limit the same result can be obtained starting from 
equation ifTTI) . 

In the Wien region of the CMB spectrum (a; ^ 1) 
again starting from H12I) and using relation Q one can 
deduce 
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(20a) 



(20b) 



Equation ((20a|l is valid for a: 1 and arbitrary (5, whereas 
(|20b(l only is applicable for x ^ 1 and x ■ 5 <^ 1. Result 
(I20b|) can be also obtained starting from equation (|17|l . 

When higher order moments are important 

Using multifrequency observations it is possible to elimi- 
nate the temperature distortion (cx {5)). In this case, the 
leading term in the signal corresponds to the y-distortion. 
Here we are interested in the question when higher or- 
der moments start to contribute significantly to the total 
spectral distortions. In the RJ region using equation l(19b|l 
one can find that for 



{5^) > I (5^) - {5^) I 



(21) 



the relative contributions of the higher moments to the 
y-distortion is less than e. This equation shows, that in 
general corrections due to the higher moments will be of 
the order of ~ | (5^) |/ (5^). For temperature distributions 
which are symmetric around the reference temperature all 
the odd moments vanish. In this case the corrections typ- 
ically will be of the order of ~ (J'*) / (^^). If we consider 
the simple case, when we compare two pure blackbodies 
with temperatures Ti and T2 and use Tret = Ti then the 
moments are given as {5^) = A*^, with A = (T2 - Ti)/Ti. 
From this we can conclude, that the relative contribu- 
tion of the third moment to the spectral distortion will 
be ^ |A|. Therefore, in the RJ region of the spectrum 
even for A ~ 10~^ higher order moments will lead to 



Figure 3. function of A for different values of 

e according to equation assuming that ((5*-') = A'^: 
For e = 10"'' (solid), e = lO'^ (dashed) and e = IQ-^ 
(dotted). The right ordinate corresponds to I'y — fcTg • 
Xy/h, with To ^ 2.725 K. 

corrections <, 1% to the y-distortion. For the CMB with 
A ~ 10~^ — 10^'^ these corrections can be safely neglected. 

In the Wien region the situation is a bit more compli- 
cated. Given some value of e one can find the frequency 
Xc above which the relative contribution of higher order 
moments will become important. To find Xc one has to 
solve the non-linear equation 



the higher multipolcs. ihcrciore we can saieiy describe 
the spectral distortions arising due to the superposition 
of blackbodies with slightly different temperatures by y- 
distortions. In what follows below we assume that higher 

- 700 order corrections are negligible. 

600 g Beam spectral distortion 

In the RJ limit the sum of blackbodies (SB) is again a 
blackbody with temperature Trj = J RTdT. The RJ 
temperature Trj can be directly measured at sufficiently 
low frequencies and hence it is convenient to compare the 
SB to a reference blackbody of temperature Trj. This 
makes any distortion due to the SB vanish at low v. 

- 400 ^ 

'°' In the following, we will call ((5 ) beam moment, when 

S = {T — Tjij)/Tjij is defined as the relative difference 
between the temperature T of a given blackbody inside the 
beam and the RJ temperature of the SB, i.e. the beam RJ 
temperature Trj. For this choice of reference blackbody 
the temperature distortion vanishes and we obtain 



52(xc)(<52) 



(22) 



or an equivalent equation directly using H12|l . In Fig.Oaic 
has been calculated for different e and A again considering 
the simple case, when we compare two pure blackbodies 
with temperatures Ti and T2 and use Trcf = Ti. Higher 
order moments contribute 1% to the total spectral distor- 
tion above Xc = 10 for A '--^ 5 • 10~^, which is roughly 5 
times the dipole amplitude. For |A| < 3 • 10~^ corrections 
are less than 1% even in the highest Planck frequency 
channels. In the limit Xc ^ 6 for estimates of the criti- 
cal frequency above which higher order correction become 
bigger than e one may use 
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(23) 



In summary, the above discussion shows that, as long 
as |A| < e the spectral distortion for x < Xc is given 
by a pure y-distortion. In this case the biggest deviations 
from a y-distortion are expected in the Wien region of 
the CMB spectrum, where the sensitivity of the Planck 
mission and other future experiments is high. But even 
for the dipole anisotropy these corrections should be less 
than 1%. If |A| ~ e deviations in the RJ become impor- 
tant and higher order corrections have to be taken into 
account. Given the amplitude of the temperature fluctu- 
ations on the CMB sky, A ^ 10^^ — 10"'^, corrections 
to the y-distortion due to higher moments will be less 
than 1% in the highest Planck frequency channels even 
for the dipole amplitude and certainly less than 0.1% for 



(24) 



where x = hv/kT^j and we introduced / = B^{T^]). We 
should mention again that the beam spectral distortion 
arising due to the SB with different temperatures can be 
approximated by a pure y-distortion with y-parameter ys 
as given by l|18|l . if only the second moment of the tem- 
perature distribution, (^^), is important. In general higher 
moments have to be taken into account leading to addi- 
tional spectral distortions. 

Sum of two Planck spectra with different weights 

For the sum of two Planck spectra with temperatures Ti 
and T2 of relative weights wi and W2 = I — wi the RJ 
temperature is given by T = wi Ti -I- W2 T2. For the first 
five beam moments one may obtain: 
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Here AT = Ti — T2 7^ was defined. Making use of equa- 
tion 1(23), in general {wi 7^ W2 7^ 0) the spectral distortion 
significantly deviates from a y-distortion by a fraction e 
at frequencies 



X > min 



3e 



12e 



1 — 2wi| V l~3wi?x;2 
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(26) 



If one of the weights wi or W2 vanishes then, as ex- 
pected, there is no distortion. In the case wi = W2, all 
the odd moments vanish and the even moments are given 
by ((5") = [AT/(Ti +T2)]". Then the critical frequency 
Xc = 2-V3 - 10-1 (Ti + T2)/|AT| above which the spec- 
tral distortion deviates significantly from a y-distortion 
becomes 10 only if |AT| « 0.069 f. This implies that for 



non gaussian temperature distributions deviations irom 
a y-distortion will only be important if the temperature 
difference becomes of the order of a few % of the RJ tem- 
perature. This is clearly not the case for the CMB in the 
observed universe. Even the dipole anisotropy leads to a 
y-distortion up to very high frequencies and we can safely 
neglect the contributions of moments of the temperature 
distribution with k > 2. 

Difference to an arbitrary reference blackbody 

Here we want to derive the equations describing the com- 
parison of the spectrum of the SB with an arbitrary refer- 
ence blackbody. This example can be applied for absolute 
measurements of the CMB. Starting again from equation 
H17() and using the results obtained for the beam spec- 
tral distortions, one can rewrite the moments of S^d = 
(T — T'rcf)/Trcf in terms of the beam moments, {S^} and 
Arj = (Trj - ri.cf)/T'i.cf, with (5 = (T - Trj)/Trj: 

- Arj (27a) 
{6^,) = (S^) + Al, + O {{6') ■ An,) . (27b) 

Using formula H17|l . one can directly infer the relative dif- 
ference between the spectrum of the SB and the chosen 
reference blackbody /rcf = B^{T^ci)'- 

= gi(iref) Arj + g^ix,^) [Al, + {5^)] , (28) 

-frcf 

where Xrcf = hv/kT^-d and we assumed that the contri- 
butions of (5^) • Ar j and any higher moments of both 5 
and Arj are negligible. Note that 52(iref) = g2{x) ■ [1 -I- 
O(Arj)] and therefore s(2(iref) (<5^) = A///, where AI/I 
is given by equation (|24|l . 

Equation H28|l shows that the spectral distortion of the 
SB spectrum with respect to the some reference blackbody 
/rcf has the following three contributions: 

(i) A temperature distortion cx Arj . This distortion van- 
ishes if the SB is compared to a blackbody with tem- 
perature Trcf = Trj. Since the y-distortion vanishes 
in the RJ region of the CMB spectrum Arj can be 
directly measured at low frequencies. 

(ii) A y-distortion cx A^^j due to the dispersion of the 
RJ temperature of the SB with respect to the cho- 
sen reference temperature. This contribution also van- 
ishes if the chosen reference blackbody has tempera- 
ture Trcf = Jrj- 

(iii) A y-distortion cx ((5^) due to the dispersion of the tem- 
perature T with respect to the RJ temperature of the 
SB. Only if the temperature distribution function is a 
Dirac (5-function, i.e. when all the blackbodies inside 
the beam have the same temperature or equivaleiitly 
when the angular resolution of the experiment is infi- 
nite, this term vanishes. For the SB it sets the minimal 
value of the y-parameter for a given temperature dis- 
tribution function. 

Equation H28|l represent the most general case for the 
comparison of the SB with any reference blackbody in the 
limit of small temperature fluctuations around the tem- 
perature of the reference blackbody. It implies that in this 



limit the y-distortion arising due to the comparison ot the 
SB spectrum with some chosen reference blackbody with 
temperature Trcf is completely determined by the second 
beam moment and the beam RJ temperature. 

3.2. Superposition of Planck spectra 

Using the results obtained so far, it is possible to gen- 
eralized the sum of Planck spectra to the superposition 
of Planck spectra. This case applies for differential mea- 
surements of the CMB, where we directly compare two 
identical beams with intensities /j and I2 and measure 
the intensity difference AI = /2 — /i. If we want to ap- 
ply the relation (O to relate AI ^ AT then we have to 
fix the reference temperature Trcf and thereby define the 
reference blackbody /rcf. During the discussion in Sect. 
12.21 we used Trcf — T (cf. equation JHJ) but in prin- 
ciple one is free to choose any temperature for which 
the limit of small spectral distortions is justified, since 
AI = I2 — Ii = I2 — /rcf + /rcf — h- In the case of the 
CMB, one will usually set Trcf equal to the whole sky 
mean temperature Tq, but as will be shown below the in- 
ferred y-distortion strongly depends on this choice. 

Defining A,; = (T^ — Trcf)/Trcf, where Ti are the tem- 
peratures of the blackbodies inside the beam i, and ap- 
plying formula (|17(l one may obtain 
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(29) 

(30a) 
(30b) 



where {X)-^^ ^ denotes the corresponding beam averages 
and we defined the beam RJ temperatures Jb,i = (^i)b i- 
In equation (|29|l we have neglected the contributions of 
any higher moments of A.^. Here the temperature distor- 
tion is directly related to the relative difference of the RJ 
temperatures. 

The y-parameter ysup,i2 can be rewritten using the 
beam moments ((5^)^ with Si — (Ti — T\^^i) /Tx^^i and 
defining A^ci,i = (A^),^ , = - T^cd/T^d- 



?/sup,12 — 2 



b,2 



(^?)b,l 



+ - [Al,^, ~ Al,^,] + 0{{5^)^^^ ■ Arcf.) . (31) 
The difference of the beam moments (5]), . does not de- 

\ ^ / b,i 

pend on the chosen reference temperature Trcf and there- 
fore always leads to the same contribution 
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(^2^)b,2 - (^?)b.l 



(32) 



to the inferred y-parameter. In the limit of high angu- 
lar resolution it vanishes. For the second term in H31|l 
(fx A^^f 2 - A?cf,i) one finds 
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(33) 



with Arj = (Jrj - Jref j/Jref and i Rj = (ib,l + ^b,2j/^- 

This contribution to the total inferred y-parameter, 
2/sup,i2(T'rcf) = 2/b + J/i2(rref), depends on the choice of 
Tref. Here it becomes clear that setting Tiof = Trj, i.e. 
to the average RJ temperature of the sum of both beams, 
only the difference of the second beam moments leads to a 
y-distortion. Let us note here, that equation H28() is a spe- 
cial case of (|29|) . where beam 1 has no contribution of the 
beam moment {Sf^^^ ^, i.e. where the beam spectrum is a 
blackbody with temperature Tb.i, and where T^^t — 7b, i. 

Dependence on the reference temperature 

Here we want to clarify some aspects of the dependence 
of the inferred y-parameter on the chosen reference tem- 
perature. If we assume that j/b <C 1 and Tb ^ > 2^,1, we 
can find the reference temperature at which the total y- 
distortion vanishes: 
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(34) 



In the limit of high angular resolution (j/b = 0) we find 
that Ty=o = Trj. One encounters this situation in dif- 
ferential measurements of the CMB brightness between 
two points on the sky. For T'b,2 > Th,! and yh > the 
reference temperature at which the y-distortion vanishes 
is always positive and also the temperature distortion 
corresponding to the first term in equation (|29|l is posi- 
tive. Then for the inferred y-distortion the following three 
regimes can be defined: (i) for T^ef < Ty=o the inferred 
y-distortion is positive, (ii) for T^ci — 7j,=o per definition 
there is no y-distortion up to first order in ?/b and (iii) 
for Tref > Ty=o the inferred y-distortion is negative even 
though 713,2 > Tb,i- This shows that for the superposition 
of Planck spectra it is in general not sufficient to choose 
the beam intensity with higher RJ temperature as I2 in 
order to obtain a positive y-parameter. 

Now we are interested in the change of Asup,i2 and 
J/sup,i2, if we change from one reference temperature Tj-ef 
to another T^f^f. For T^.^^ one can again write down an 
equation similar to 1)29(1 . where now = (Ti — T^ei)/Tfci 
is replaced by = (T, - T4f)/r4f and x = hv/kT.^i by 
x' — hi'/kT'.^f. Using (|30|) and H31|) . we can write 
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^sup,12 ' ^RJ 
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with Aj- — [T-rcf — T^ci) /'^ici ^^^^ where prime denotes, that 
the reference temperature T/^^f was used. The difference of 
y-parameters is given by Ay = y^up^iz - ysup,i2 = Asup,i2 • 
A].. For Ty-ci = Trj the contribution yi2 = 0. Therfore 
we can write the inferred y-parameter for any reference 
temperature as: 

X , Jb,2 - Jb,l Tb,l + Tb,2 - 2Trof , , 

y(Trcf) « 2/b + — . (36) 



Th,l + lb, 2 



ihis result describes all the cases discussed previously: It 
we want to discuss the beam distortions we have to set the 
reference temperature to Tr j making ?/ « yb- If we are dis- 
cussing the y-parameter for the comparison of the SB with 
some reference blackbody we have Tb,i = Tiof resulting in 
y ~ J/b + A|,j/2. This implies that in the limit of small 
spectral distortions the y-parameter for the superposition 
of Planck spectra is completely determined by the second 
beam moments and the beam RJ temperatures. 



4. Superposition of two Planck spectra 

In order to illustrate the main results obtained in the pre- 
vious Sect, we now discuss the spectral distortions arising 
due to the sum and difference of two pure blackbodies with 
different temperatures Ti and T2 and the corresponding 
intensities Ii — B^{Ti) and I2 — B^[T2) in the presence 
of a reference blackbody /icf = B^{T-[ci) with tempera- 
ture Ticf. In principle one can directly apply the results 
obtained in Sect.|31and easily derive the equations describ- 
ing this situation, but in order to check the derivations of 
the previous Sect, we here choose an approach starting 
from the expansion taking only first order corrections 
into account. Since the CMB temperature fluctuations are 
very small this approximation is sufficient. 

The results obtained for the sum of two blackbod- 
ies apply to the case of an absolute measurement of the 
CMB sky, where the beam contains only blackbodies with 
two equally weighted temperatures. In the limit of narrow 
beams the results obtained for the difference of two black- 
bodies can be directly used to discuss the effects arising 
in differential measurements. In Sect. IHlwe will use some 
of the result of this Sect, in the discussion about cross 
calibration issues. 



4.1. Sum of two Planck spectra 

We want to describe the difference A/s — Is — lid, as 
a function of frequency. Here we defined the sum of the 
two blackbodies by /s = {Ii + l2)/2. The easiest is to 
first calculate the difference A/j = Ii — I^ci since A/s = 
(A/i + A/2)/2. If we define A, = {Ti - r„f)/rrcf and 
X — hv/kTid we can make use of equation Q and write 
the relative difference of /s and I^cf as 
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Here we abbreviated the y-parameters as j/i = A^/2. 
Defining the RJ temperature Trj = (Ti + r2)/2 of the 
sum of the two blackbodies, the sum of the first order 
terms and the y-parameters may be rewritten like 
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From this it can be seen that if the reference temperature 
is set to Trj the temperature distortion corresponding to 




Figure 4. Dependence of the inferred y-parameter on the 
chosen reference temperature Tj-^f for two blackbodies with 
temperatures Ti = 2.72162 K and T2 = 2.72838 K, cor- 
responding to the minimum and maximum of the CMB 
dipole. In this case the mean RJ temperature is Trj — 
(Ti + T2)/2 = 2.725 K and yb = 7.7 • lO"'^. — Sum of two 
blackbodies: total y-parameter as given by (|38b|) (dashed 
line). — Difference of two blackbodies: the absolute value 
|?/dif I of the inferred y-parameter as given in equation (|41bp 
(solid line). The dotted vertical lines show the location of 
Trof = Ti and Ticf = T2 respectively. 

the first term in (|37() vanishes (Arj — 0) and only a min- 
imal y-distortion with y-parameter 



T1-T2 



T1+T2 



(39) 



is left. At low frequencies (x < Xy) this reflects the fact 
that the sum of RJ spectra is again a RJ spectrum. 
Comparing (|38|l to (|27|l one may conclude that j/b cor- 
responds to the contribution of the second beam moment 
((5^) to the y-parameter. Comparing (|39|l to (|25b(l with 
wi = W2 = confirms this conclusion. 

The dependence of the inferred y-parameter on the 
chosen reference temperature is illustrated in Figure 0] 
As mentioned above the y-parameter reaches a minimum 
for Tict = Trj. It is important to note that the min- 
imal value of the y-parameter j/b does not depend on 
the reference temperature Trcf but only on the second 
beam moment of the temperature distribution function. 
Changing Ticf from Trj — > Ti the inferred y-paramcter 
monotonically increases being a factor of ~ 2 larger than 
?/b for Trcf — Ti. This can be understood as follows: If 
we set Tref = Ti we obtain ysum = [Ti - T2Y / ATf k. 
(Ti - T2f/AT^ + 0{/\T^/T^), with AT = Ti - T2. 
Therefore the ratio y.um/yd « (^i -f T2f /2T^ w 2. 

Outside the region Ti < Ticf < T2 the inferred y- 
parameter increases further. It is in principle possible to 
gain large factors in the inferred y-parameter by going very 
far away from the RJ temperature of the SB, but at some 
point next order corrections will become important and 
therefore limit T^cf to the region where the approximation 
of small spectral distortions is still valid. For the CMB 
the most natural choices of the reference temperature are 



the whole sky mean temperature Jq a-nd the maximum 
or minimum of the CMB dipole. For the SB the inferred 
y-parameter varies from y ^ 7.7 • 10^^ for Trcf = Tq to 
y - 1.5 • 10"*^ for Tref = Ti. The behavior of the inferred 
y-parameter shows, that in order to minimize the arising 
spectral distortion for absolute measurements it is impor- 
tant to use an internal calibrator with temperature close 
to the beam RJ temperature Trj. 

In the limit Ti — > T2 we are in principle comparing one 
pure Planck spectrum with a reference blackbody. This 
case applies for an absolute measurement of the CMB sky 
in the limit of very narrow beams. In this case it follows 
that yb and in addition ysum ~> for Tj-of = Trj. 
Except for some changes in scales the behavior of the curve 
shown in Fig. 01 is unaffected. 

4.2. Difference of two Planck spectra 

Since AI = I2 ~ h = A/2 - A/i, where A/; = 1; - I^ci, 
applying equation one can easily find 
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Again it is possible to simplify the difference of the first 
order terms and the y-parameters leading to 
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where Arj = (Trj -T„f)/Tref and Trj = (ri+T2)/2. For 
the difference of two Plancks the temperature distortion 
does not vanish for any choice of the reference temperature 
(see equation H41a(l ^. whereas the y-distortion changes sign 
at Trcf — Trj. As mentioned earlier, this shows that for the 
difference of two Planck spectra in order to obtain a posi- 
tive y-parameter it is in general not sufficient to choose the 
intensity with higher RJ temperature as 12- Comparing 
(|41|l to (|3()a|l and 133|l and keeping in mind that both beam 
moments {Sf^^^ ^ = shows the equivalence of both ap- 
proaches. In Figure^the dependence of the absolute value 
of the inferred y-parameter on the chosen reference tem- 
perature is shown. At Ticf = Ti the inferred y-parameter 
for the difference of two blackbodies is 2 times larger than 
in the case of the sum, (|j/dif |/2/sum)| j-. ~ 2. 

As in the case of the sum of two Plancks for Tjof out- 
side the region Ti < T^ef < T2 the inferred y-parameter 
increases strongly. Again there is a limit set by the approx- 
imations of small spectral distortions. The behavior of the 
inferred y-parameter shows, that in order to minimize the 
spectral distortion arising in differential measurements it 
is important to choose a reference temperature close to 
the RJ temperature of the observed region. 



5. Spectral distortions due to the CMB dipole 

The largest temperature fluctuation on the observed 
CMB sky is connected with the CMB dipole. Its am- 
litude has been accura tely measured by Cobe/Firas 
iFixsen fc Matheill2002 l): ATd = 3.381 ± 0.007mK, corre- 
sponding to AT/T ~ 0.2 % on very large angular scales. 



Assuming that it is only due to the motion oi the so- 
lar system with respect to the CMB rest frame it im- 
plies a velocity oi v — 372 ± 1 km/s in the direction 
{l,b) = (264.14° ± 0.15°, 48.26° ± 0.15°). In order to un- 
derstand the spectral distortions arising due to the dipole, 
we start with the direction dependent temperature of the 
CMB, neglecting any intrinsic anisotropy: 

To 



7[1-/3m] 



i 



To 



1 - 



6 



(42a) 



(42b) 



Here we introduced the abbreviation /i = cos 9, with 9 
being the angle between the direction of v and the location 
on the sky. Tq denotes the intrinsic CMB temperature, 
P — v/c is the velocity in units of the speed of light and 
7 = — P'^ is the corresponding Lorentz factor. The 

whole sky mean RJ temperature is given by 
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where 'f denotes full sky average temperature. Equation 
1)43(1 implies that due to the motion of the solar system rel- 
ative to the CMB restframe the observed whole sky mean 
RJ temperature is AT = 0.70 /iK less than the intrin- 
sic temperature Tq. The measured value of the whole sky 
mean temperature as given by the Cobe/Firas experi- 
ment l|Fixsen fc M ather. 2002) is To = 2.725 ± 0.001 K. 
Therefore this difference is still ~ 1000 times smaller than 
the current error bars and far from being measured. 

Below we now discuss the spectral distortions arising 
due to the CMB dipole in the context of finite angular 
resolution. Since only the second moment of the temper- 
ature distribution is important, we will use the second 
order expansion of T{fi) in /3 as given in equation 142b(l . 
In the last part of this section we show that, starting with 
equation (|42a|l all the results of this section can also be 
directly obtained by expansion of the blackbody B^{T{^)) 
in terms of small (3. This shows the equivalence of both ap- 
proaches. Nevertheless, the big advantage of the treatment 
developed in Sect. |3 is that it can be applied to general 
temperature distributions and that the source of the y- 
distortion can be directly related to the second moment 
of the temperature distribution. 

5.1. Whole sky beam spectral distortion 

Defining 5 = [T(yLt) — TfJ/Tf, where Tf is given by equa- 
tion the full sky moments can be calculated by the 
integrals {S'')^ = 5 /^-^ S{ii)'' d/i. The first three moments 
may be found as 

{S°)^ = l, {S'), = (44a) 
('5')f = ^-i'J y + <^(/3')- (44b) 

where 5 = \n[j^] . With /3d = ATd/To = 1.241 • lO^^ one 
finds: ((5^)f — 5.1-10^^. This implies a full sky y-distortion 
with 



yd = f «2.6•10-^ 
6 



(45) 



which IS currently ~ 60 times below the COBE/i" IRAS up- 
per limit. This distortion translates into a temperature 
difference of AT = 0.7 • gy{x) /xK. The exact behavior of 
gy{x) is shown in Fig. [3 We have checked that devia- 
tions from a y-distortion will become important only at 
extremely high frequencies {x > 360). 

Here we want to note, that this full sky distortions 
arises due to the existence of the dipole anisotropy. In 
absolute measurements of the CMB ?/d places a lower limit 
on the full sky y-parameter (see Sect. I3.ljl . 

5.2. Beam spectral distortion due to the CMB dipole 

Here we are interested in the angular pattern of the y- 
distortion induced by the SB over the dipole for an obser- 
vation with a finite angular resolution and in particular in 
the location of the maximal y-distortion, when we com- 
pare the beam flux at different frequencies to a reference 
blackbody with beam RJ temperature. To model the beam 
we use a simple top-hat filter function 
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where the z'-axis is along the beam direction^ and 9^ is 
the radius of the top-hat in spherical coordinates. 

The fcth moment of some variable X over the beam is 
then given by the integral 
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where 'r' denotes beam average for a top-hat of radius 
and /ir — cos 9,- . Defining fio — cos 9o and using 
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one can calculate the beam averages of the dipole and 
quadrupole anisotropy 
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where we introduced the abbreviation (± ~ These 
equation will be very useful in all the following discussion. 

Using formula l(42b|l and equations 1)47(1 and ((49(1 the 
mean RJ temperature inside the beam in some direction 
(00 J ^0) relative to the dipole axis is given as 



T, = To 



1 - Y + /3C+ Mo + /9Vr C+ ( Mo - 3 



(50) 



The position of the maximum and minimum is, as ex- 
pected, at 00 = and 9o = respectively. 

Now the second moment oi 5 = [T(/i) — Ti-J/Ti. inside 
the beam can be derived as: 
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^ Prime denotes coordinates with respect to the system S' 
where the z'-axis is defined by the direction of the beam. 
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Figures. Angular distribution of j/r as given by equa- "-"^ ^ / 
tion H52|l for 9i- — 20°. The maximal spectral distortion is _ / 

expected to appear in a ring perpendicular to the dipole / 

axis and in this case has a value j/r.max = 2.3 • 10^^. The "'^ ^ / 
minimal y-parameter has a value of y = 2.4 • 10^^'' and is - / 

located around the maximum and minimum of the dipole. y 

Q 1 — . ^. 1 r' — \ 1 1 1 1 \ 1 1 1 1 \ — 

The galactic plane was cut out using the kpO-mask of the o so loo iso 

Wmap data base. 



where up to second order in /3 only the CMB dipole con- 
tributes and we have set /3 = /3d. The corresponding y- 
parameter 



2/r = yd 
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(52) 



has a monopole and quadrupole angular dependence. As 
expected it vanishes for high angular resolution (^_ — > 0). 
In Fig.Elthe angular distribution of j/r for 6^ = 20° is illus- 
trated. The maximum lies in a broad ring perpendicular 
to the dipole axis, which even without taking the galaxy 
into account is covering a large fraction of the sky. In the 
shown case, the maximal value of the y-parameter is ap- 
proximately 2/r, max ~ 2/d/lO. 

Using equation H51|l the position of the maximum can 
be found: 0o,max — f ■ This suggests that the location 
of the maximal distortion is where the derivative of the 
temperature distribution is extremal. The maximal y- 
parameter is given as 
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The dependence of j/r.max on the beam radius is shown 
in Fig. ini For beam radii larger than 9^ = 90° it follows 
that ?/r,max > Vd and at O^- = 120° it is 9/8 times yd- The 
maximal y-parameter, ?/r,max, is a steep function of beam 
radius: It has values yd/2, yd/3 and yd/10 for 9r — 51.83°, 
9r = 40.68° and 9^ — 21.29° respectively and vanishes like 
~ 3 0^/4 for high angular resolution. For = 5° it is ~ 6 • 
10"'^ X ?/d- This implies that the beam spectral distortions 
resulting from the dipole are lower than AT/T ^ 10~^ for 
experiments with resolution better than ~ 10°. As will be 
shown later, on this level contributions from the higher 
multipoles will start to become important (Sect. El). 



Figure 6. Dependence of the y-parameter in units of yd = 
2.6 • 10^^ on the beam radius 9^ in degree for a circular 
beam: for /io = ±1 as given by equation H52|l (dash- 
dot-dotted), yr,max accordiug to equation (|^ (solid). For 
values — 1 < /io < 1 the curves lie between these two 
extremes. 



second beam moment (5^) and another due to the dis- 
persion Ap J of the beam RJ temperature relative to the 
temperature of the chosen reference blackbody (cf. equa- 
tion Ij27b|l ). Therefore one can immediately write down 
the relative difference of the beam intensity and the in- 
tensity of the chosen reference blackbody making use of 
(|?7jl . f^ . ISU and In the hmit 9^ ^ we can apply 
the results obtained in Sect. l4.ll for the sum of two black- 
bodies, if we set Ti = T2. Here we want to discuss the case 
Tj-of — Tq in more detail. For this choice the residual full 
sky y-distortion (which may be obtained by averaging the 
intensity data over the whole sky) is minimal. 



Distortion with respect to Tq 

Now, instead of comparing the spectrum of the SB to a 
blackbody with beam RJ temperature T,- we compare to a 
blackbody of temperature Tq. Then the spectral distortion 
follows from equation 1)28(1 and in addition to the beam 
moment ((51(1 the first two moments of A = (Ti- — To)/To 
have to be calculated: 



Ai = 
A2 = 



/32 /I 

- Y + /3d C+ Mo + Pi A*r C+ - 3 



(54a) 
(54b) 



5.3. Distortion with respect to any Tref 

We have seen in Sect. 13.11 that the inferred y-parameter 
for the comparison of the SB with some reference black- 
body with temperature Tpef has one contribution from the 



The spectral distortion arising from the SB can now 
be characterized as follows: (i) At low frequencies the 
y-distortion vanishes. The motion of the solar system 
with respect to the CMB restframe induces a temper- 
ature dipole (~ /3d) and a temperature monopole and 



quadrupoie (~ p^j all resulting irom the A -term: 
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Pd Air 



(55) 



In the hmit of high angular resolution (^+,/Xr —^ 1) this 
is equivalent to the expansion l|42b(l . (ii) A y-distortion is 
induced which is proportional to the sum + (S^)^.'- 



AT 
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(56) 



It has a monopole and quadrupoie angular dependence 
and only arises due to the CMB dipole. The sum of both 
contributions mentioned above can be rewritten as 



AT 



(5^)^ + ^ [A^ + (5^)J (57a) 



2 



C+ Mo 



(57b) 



where we used the definition ^ for 5(2;). Rewriting the 
function (72 (i) in terms of g{x) lead to exact cancellation of 
the motion induced temperature quadrupoie and changes 
the temperature monopole by /3^/3 (cf. A — A^ — ((5^)^.- 
tcrm in equation l|57a|) '). 

Here we want to note, that in the picture of the SB it 
is easily understandable, that there is no difference in the 
resulting spectral distortion whether the dipole anisotropy 
is intrinsi c or due to m otion. Therefore, as has been noted 
earlier bv iKamionkowski fc Knox (20Q3i), it is impossible 
to distinguish the intrinsic dipole from a motion induced 
dipole by measurement of the frequency dependent tem- 
perature quadrupoie. 

Expansion of AI/Iq for small (3 

Here we want to show, that the results (|51|l and (|57|l can 
also be directly obtained starting with the expansion of 
the blackbody spectrum with temperature T(/i) for small 
velocity (3 and thereby prove the equivalence and correct- 
ness of both approaches. For this, inserting H42al) into the 
blackbody spectrum JHjl, by Taylor expansion up to second 
order of (3 one may find 



^-/3'5i(i) 
Jo 



+ /3.gi(i)/i 



(58) 



This e quation ha s be en obtained earlier by 

ISazonov fc Sun vaevI for a cluster of galaxies 

mo ving with respect to the CMB and was later applied 
by IKamionkowski fc Knoxl (|200^ to discuss aspects of 
the observed CMB dipole and quadrupoie. It is valid for a 
measurement of the CMB temperature with high angular 
resolution, where the y-parameter due to the second 
beam moment is negligible (y ^ 10~^ for 10° angular 
resolution). 

First we want to calculate the SB in a circular beam 
and compare it to the blackbody /q of temperature Tq, 



I.e. we want to derive Using equation i|47(l and 

since only /i and /i^ depend on /i' and 0' this integration 
with H49|l immediately leads to equation (|57|l . 

Next we want to derive the spectral distortion relative 
to the reference blackbody with beam RJ temperature 
starting with 158|) . For this we rearrange (|58|l in terms of 
gi{x) and 52(2;) leading to 

I{x)^Io{x)-[l+gi{x)-So+92 (x) ■ S^] , (59) 



with (5o = — 7^ 



-g- , f3 fi + 13"^ (fi'^ - I) . In the RJ limit 
gi — > 1 and .92 ^ 0. This immediately leads from equation 
dSni to Trj = (To • [1 + 6a])^. = To • [1 + A], where A is 
given by equation (|54a|l . Now we replace x = hv/kT^ — 
{hv/kT^j) ■ T^j/Tq = • [1 A] in equation and 
expand in second order of A obtaining: 

I(x^) = h(x^) ■ [1 -I- gi(x^) • ((5o - A + A^ - 5^A) 

+ 32 (x*) -(50+^" -2(50 A)] , (60) 

with /, (x*) ~ B^{T^j). Since only Sq depends on /i' and 
(/)' we may perform the beam averages using equations 
(|47|l and (|49|l . This is equivalent to replacing Sq ^ A and 
Sq ('5o)r i'^ equation (|60|l . Therefore this leads to 



/-/, 



g2(x.)-[{6l)^-A' 



(61) 



Inserting equation (|54b(l and 



(62) 



one may obtain (Sq)^ — A^ = {S"^)^, where (5^)^, is given 
by equation (|51|) and wc made use of the identities 



i-C+ = [i + C+]C- 



(63a) 
(63b) 



Let us note here, that since only Sq depends on fi' and (j)' 
we may have calculated the beam average already from 
equation (|59|l . After expansing in terms of A this would 
have directly lead to the result (|61|1 without the intermedi- 
ate step (|60|l . This finally shows the complete equivalence 
of both approaches. 



6. Spectral distortions due to higher multipoles 

The largest temperature anisotropy on the CMB sky is 
connected with the CMB dipole which was discussed in 
detail above (Sect. EJ. In what follows here we are only 
concerned with spectral distortions arising from multi- 
poles with / > 2. As has been shown above, the spectral 
distortions induced by the CMB dipole are indistinguish- 
able from a y-distortion. Since the typical amplitude of 
the temperature fluctuations for /> 2 is a factor of ^ 100 
less than the dipole anisotropy (AT/T ~ 10~^), the spec- 
tral distortion arising from higher multipoles will also be 
indistinguishable from a y-distortion. 



Lalculations using a realization ot the livid sky 

In order to investigate the spectral distortions arising from 
the multipoles / > 2, we use one realization of the CMB 
sky computed with the Synfast code of the Healpix 
distribution^ given the theoretical Ci's for the tempera- 
ture anisotropics compute d with the Cmbfast code'^ for 
the WMAP best fit model ijBennett et alll2003|) . Here we 
are only interested in the distortions arising from the pri- 
mordial CMB anisotropics. Future experiments will have 
~ arcmin resolution. Therefore we perform our analysis 
in the domain, I < 3000, corresponding to angular scales 
with dimensions 6 ^ 3.6'. 

Using the generated CMB maps we extract the tem- 
perature distribution function R(T) inside a circular beam 
of given radius 9^ in different directions on the sky. We 
then calculate the spectrum and the spectral distortion 
inside the beam using equation 1)12(1 and setting the ref- 
erence temperature to the beam RJ temperature Trj = 
/ RTdT. As expected, we found that the spectral dis- 
tortions are indistinguishable from a y-distortion with y- 
parameter given by equation 1(181) . As discussed above this 
follows from the fact, that the temperature fluctuations for 
I > 2 are extremely small (~ 10"^). 

In Fig.[3the probability density p{y) for a given aper- 
ture is shown. It is defined such that p{y) dy gives the 
probability to find a value of the y-parameter in a random 
direction on the sky between y and y + dy for a given beam 
radius. As expected, the average value (y) increases with 
aperture radius, while the width of p{y) decreases. For 
9i ^ 1 the average beam y-parameter varies only slowly. 
In the limit of a whole sky measurement, the y-parameter 
converges to yh = 8.33 • 10^^", which corresponds to the 
whole sky rms dispersion AT ~ 111 /iK of the used real- 
ization of the CMB sky. For O^- <^ 1 the beam y-parameter 
falls off very fast. As Fig. 13 suggests, for 9r ^ 2 — 5° the 
distortions arising due to the dipole arc negligible in com- 
parison to the distortions arising due to higher multipoles. 

FigureHalso shows the y-parameter resulting from the 
dispersion of the beam RJ temperatures in different direc- 
tions on the sky relative to Tq. It drops with increasing 
beam radius, because the average beam RJ temperature 
becomes closer to Tq for bigger beam radius. Looking at 
the sum of the beam distortion and the distortion due to 
the dispersion of the beam RJ temperatures shows that 
the averaged whole sky distortion is independent of the 
angular resolution. 

In Fig.|Slthe cumulative probability for finding a spec- 
tral distortion with y > ya is shown for different apertures. 
Fixing the value yo, there is a maximum probability of 
2.4% to measure y > 10^^ for an aperture radius 9^ ^ 1° 
during the mapping or scanning of extended regions of the 
sky. This correspond to ^ 500 sources on the whole sky. 
In Fig. 13 the dependence of the number of sources with y 
above yo for a given radius 9^ is illustrated. It peaks around 
9i- ~ 0.5 degree, corresponding to the first acoustic peak. 
We should mention that the maxima of the y-distortion 
do not coincide with the maxima of AT/T, but that they 
both should have similar statistical properties. We expect 

^ http:/ /www. eso.org/science/healpix/ 
^ http://www.cmbfast.org/ 




Figure 7. Top: Probability density p{y) for different aper- 
ture radii: 9, = 0.25° (sohd), 9, = 0.5° (dashed) and 
9i- — 1.0° (dashed-dotted). — Bottom: Dependence of 
the average y-parameter on the aperture radius 9,- in de- 
grees: for the beam spectral distortion (solid) and for the 
y-parameter resulting from the dispersion of the beam RJ 
temperatures in different directions on the sky relative to 
To (dashed). Also shown is the maximal beam spectral 
distortion arising due to the dipole, ?/r,max (dotted), ac- 
cording to equation . 



that the maxima of the y-parameter are there, where the 
derivatives of the averaged temperature field are large. 

The motion induced CMB octopole 

In section [S] we discussed spectral distortions due to the 
motion of the solar system relative to the CMB restframe 
up to second order in (3^- On the level of y ~ 10"^" — 10~^ 
third order corrections cx (3^ should also start to con- 
tribute. In this order of /3d not only the motion induced 
dipole and quadrupole lead to spectral distortions but 
also terms connected to the products of the intrinsic 
dipole and quadrupole with the motion induced dipole and 
quadrupole. Here we now only discuss the motion induced 
CMB third order terms. 

Expanding ((42a(l up to third order in /3d and calcu- 
lating the beam averaged RJ temperature for the top-hat 
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Figure 8. Top: Cumulative probability p{y > yo) for dif- 
ferent aperture radii 9,- = 0.1° (solid), 9^ = 0.25° (dashed) 
and 0r = 1.0° (dashed-dotted). - Bottom: Dependence of 
the cumulative probability p{y > ya) on the aperture 9^ 
for 2/0=8- 10-1° (sohd), yo = 9 ■ lO^^" (dashed) and 
yo = 10^^ (dashed-dotted). 
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Figure 9. Number of regions on the sky with y > yo for 
different aperture radii 9^ in degree: for yo = 8 ■ 10"^° 
(sohd), 2/0 = 9- 10^1" (dashed) and yo = 10"^ (dashed- 
dotted). 



moment to the y-parameter. If we consider the third beam 
moment we find 



(5^)^ = 6 yd /3d 
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C+C- (66) 



It reaches a maximal value of {S^}^ /6 ~ 1.8 - 10"^^. Due 
the strong dependence on 9r we may therefore also neglect 
any contribution of the third beam moment. 



beam using 147|l leads to 



'7r,oct —Tr + To (3^ C+ 
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(64) 



where Ty is given by 1)50(1 . In addition to the motion in- 
duced octopole in third order of /3d there also is a correc- 
tion to the CMB dipole. Now we can also recalculate the 
second beam moment leading to a y-parameter of 



yr,oct = ydPd- 



Mo 



^Mr^ - e^io^Mo - 



C+C- 



(65) 



for the third order in /3d- The total motion induced y- 
parameter is then given by the sum y = yr + yr,oct, 
where y^ is given by (|52() . As before in the limit 6'r — > 
the y-parameter yr,oct vanishes. It has a dipole and oc- 
topole angular dependence and reaches a maximal value of 
yr.oct 9.6 - lO""'^^ in the /io, 6'r-plane. Therefore one may 
completely neglect the contribution of the second beam 



7. Spectral distortions induced in differential 
measurements of the CMB sky 

In Sect. we discussed the spectral distortion arising due 
to the CMB dipole inside a single circular beam in com- 
parison to a reference blackbody with temperature Tq. In 
this chapter we address the spectral distortion arising in 
differential measurements of the CMB fluctuations, where 
two beam intensities Ii and I2 are directly compared with 
each other and the intensity difference A/ = I2—I1 is mea- 
sured. Since in Sect.lHlwc have shown that y-distortion 
arising due to higher multipoles have corresponding y- 
parameters y < 10~^, here we are only taking distortions 
arising due to the CMB dipole into account. If we assume 
that the both beams are circular and have the same ra- 
dius 9r (see Fig. irn]) . we may define the beam RJ temper- 
atures as Tr^i and where T^^i of each beam is given 
by equation ifsUI) . Now, defining the reference blackbody 
/* = By(T^) with temperature T^ef = 7* using equations 
and 1(29(1 we find for the inferred temperature difference 
at frequency i* — hu/kT^ei 



AT(x^) 



Ai.,12 -I- yr,12 - ffa(i*) , 



(67) 




Figure 10. Illustration of a differential observing strat- 
egy: The maximum of the dipole lies in the direction of 
the z-axis. Both observed patches have the same radius 6^ 
and the observed intensity difference is A/ = I2 — Ii- 



where we have introduced the abbreviations 



Tr,2 — Jr.l 



2/r,12 = yb,r + Vrcf 

Ti-2 — Tr,i Ti-j + Tr.2 — 2 T* 



Vrcf 



(68a) 
(68b) 
(68c) 



Using equation l|52|) . the difference of the beam y- 
parameters yb,r = J/r(/^2) — yr(Mi) can be written as 



-3 yd A/i [2 + A/z] C+ C- , 



(69) 



where /i^ = cos 9i, A/i = /i2 — /^i and 6i is the angle be- 
tween the dipole axis and the beam i. Here it is important 
to note, that in the case A/i = 2, i.e. when we are com- 
paring the maximum and minimum of the CMB dipole 
2/b,r = 0, since both beams have the same shape and size. 
In this case we are directly dealing with the difference of 
two blackbodies with temperatures Tr.i and we may apply 
the results obtained earlier in Sect.0] 

Now one can write for the relative temperature differ- 
ence of the two beams at two frequency and x^^2, with 
a;*^2 > and x^,^i = hvi/kT^, as 



AA 



Ar(x,,2)- Ar(x,,i) 



(70) 



If we assume that x■^,^l <C 1, i.e. for a measurement in the 
RJ region of the CMB spectrum, then we can approxi- 
mate 3y(a;*,i) « x1 i/Q. For estimates we will assume that 
gy{x*,i) = 0. 

On the CMB sky the most natural choices for the ref- 
erence temperature is the full sky mean temperature Tq 
and the temperature of the maximum or the minimum 
of the CMB dipole anisotropy, whereas the last two are 
equivalent. Below we now discuss the dependence of the 
inferred y-parameter on the angle between the beams for 
these two choices of the reference temperature. In the limit 
0r — ^ we can apply the results obtained in Sect. 14.21 for 
the difference of two blackbodies. 
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Figure 11. Dependence of the y-parameter in units of 
yd = 2.6 • 10~^ on the beam radius 9^ in degree for dif- 
ferential measurements: y^.o.opt as given by equation H73|l 
(dashed/left ordinate), y^,i,opt as given by equation H76|l 
(solid/right ordinate) and yA,2,opt as given by equation 
(TTTjl (dotted/left ordinate) . 

Case = Tq 

In this case, making use of equation (|5()|l and H68c|l one 
may find 

2/,ef,o - 3 yd A/i [2/ii + A/i] d+Oifil). (71) 

Using the beam y-parameter yb.r as given in equation H69|l 
we find for the corresponding total inferred y-parameter 



VAfi - 3 2/d iMr C+ A// • [2 /ii + A/i] I + 0{l3l) 



(72) 



The largest values of the y-parameter in the /ii , A/i- 
plane are obtained for the combinations (/^i, A/i) = (0, 1), 
(0,-1), (1,-1), (—1, 1). In these optimal cases it follows 



2;^,o,opt = 3ydC+ l/^r| +0{Pl) . 



(73) 



The dependence of y^.o,opt on the radius of the beam is 
shown in Fig. ^] The distortion does not vanish for high 
angular resolution {O-c — > 0) and the maximum is 3 times 
bigger than yd- The corresponding absolute temperature 
difference in some of the Planck frequency channels are 
given in Tabled For 9,- = 7r/2 the y-distortion vanishes. 

Let us note here, that in the case A/i = 2, i.e. when 
we are directly comparing the maximum and the mini- 
mum of the CMB dipole, the spectral distortion vanishes. 
This can be understood as follows: As was argued in Sect. 
13.21 if the chosen reference temperature is equal to Ty^Q 
as given by equation H34|l the total y-parameter is zero. 
For /ii = —1 and /i2 = 1 the beam moments {Sif)^ ^ 
are both equal and therefore do not contribute to the 
total inferred y-parameter Ij68b|) . In this case we obtain 
Ty=o = Trj = (T.a-f T;^2)/2 = To+O{l3l). Therefore, con- 
sistent in second order of (3d the total inferred y-parameter 
vanishes. This conclusion can also be directly drawn from 
equation 



Lase 1^ = ir,i 

In this case, again making use of equation (|5()|l and H68c|l 
one may find 

2;ref,i = 3ydA/i2^2 _^0(/33). (74) 

Now using equation (j69|l we may write the total inferred 
y-parameter as 

yA,i = 3ydC+ I Am • [/^rM2 - Mill + 0{PI) . (75) 
The largest values of the y-parameter in the /ii, A/i-plane 
are obtained for the combinations (/ii,A/i) — (1,2) and 
(— 1, —2). In these optimal cases the y-parameter is 

2/Ai,opt = 12ydC+ + 0(/3^1). (76) 

For high resolution (C+ — > 1) we can obtain y^^i.opt di- 
rectly as the difference of two Planck spectra with tem- 
perature corresponding to the maximum T2 and the min- 
imum Ti of the dipole using (|41b(l with Tj-of — Ti. In this 
case we directly get \ [{T2 - Ti)/Tif = 2(31 = l^yd- As 
will be discussed later this way of comparing the max- 
imum and minimum of the CMB dipole should open a 
way to cross calibrate the frequency channels in future 
experiments with full sky coverage to very high precision 
(Sect. El). 

The dependence of j/^,i,opt on the beam radius is 
shown in Fig. |S| As in the previous case, the y-distortion 
does not vanish for high resolution {0^ 0) and in the 
maximum it is even 12 times bigger than j/d, correspond- 
ing to j/opt = 3.1 • 10^^. This is only ~ 5 times below 
the current upper limit on the mean y-parameter given 
by Cobe/Firas. The corresponding absolute tempera- 
ture difference in some of the Planck frequency channels 
are given in Tabled 

Another case may be interesting for CMB missions 
with partial sky coverage, but which have access to the 
maximum or the minimum of the dipole (/^i = ±1) and 
a region in the directions perpendicular to the dipole axis 
{^i2 = 0): Using (|75|) with (^i, A/i) = (±1,0) one obtains 

yA2,opt = 3ydC+ + 0(/3|). (77) 

In this case the amplitude is comparable to 2/^,0 1 but the 
dependence on the beam radius O-c is much weaker (see 
Fig.EJ. 

8. Cross Calibration of frequency channels 

The cross calibration of the different frequency channels 
of CMB experiments is crucial for the detection of any fre- 
quency dependent signal. As has been mentioned earlier 
usually the dipole and its annual modulation due to the 
motion of the earth around the sun is used for calibra- 
tion issues. The amplitude of the dipole is known with a 
precision of ^ 0.3% on a level of AT ~ niK. But the sen- 
sitivities of future experiments will be significantly higher 
than previous missions and request a possibility for cross 
calibration down to the level of tens of nK. 

In this Sect, we now discuss two alternative ways to 
cross calibrate the frequency channels of future CMB ex- 
periments. Most importantly we show that making use of 
the superposition of blackbodies and the spectral distor- 
tions induced by the CMB dipole should open a way to 
cross calibrate down to the level nece ssary to detect sig- 
nals from the dark ages as proposed bv lBasu et al.. (,2003,) . 



b.l. Lalibration using clusters ot galaxies 

It is obvious, that brightest relaxed clusters of galax- 
ies on the sky open additional way to cross calibrate 
the different frequency channels of future CMB experi- 
ments like Planck and ground based experiments like 
Act. The CMB signal of the majority of rich clusters has 
very distinct y-type spectrum due to the thermal SZE. 
Multifrequency measurements allow to determine their y- 
parameters and thereby open a way to use them as sources 
for cross calibration. This cross calibration will then en- 
able us to detect the effects of relativistic correction for 
very hot clusters and thereby to measure the temperatures 
of clusters independent of X-ray observations. 

Clusters such as Coma, where existing X-ray data pro- 
vide excellent measurements of the electron densities and 
temperatures inside the cluster and where we know the 
contamination of the CMB brightness by radiosources and 
the dust in galaxies very well, even may be used for abso- 
lute calibration of CMB experiments. For the Coma clus- 
ter today X-ray data allows us to predict the CMB surface 
brightness with precision of the order of a few percent. 

Due to the redshift independence of the SZ sig- 
nal, distant clusters have similar surface brightness but 
much smaller angular diameters. Therefore they are good 
sources for calibration issues for experiments covering only 
limited parts of the sky and having high angular resolu- 
tion like Act (for example A ct is p laning to investigate 
200 square degree of the sky ijKosows kv. 2003)). 

One disadvantage of clusters for cross calibration pur- 
poses is that they are too bright: Typically clusters have 
y-parameter of the order of ^ 10^^ — 10"'' and therefore 
are only a few times weaker than the temperature signal 
of the CMB dipole 3mK). In comparison to the level 
of calibration and cross calibration achieved directly using 
the dipole this is not much of an improvement. 

In addition planets or galactic and extragalactic ra- 
diosources might serve for cross calibration purposes, if 
we know their spectral features with very high precision. 
Quasars and Active Galactic Nuclei are not good candi- 
dates for this issue, since they are highly variable in the 
spectral band of interest. 

8.2. Cross calibration using the superposition of 
blackbodies 

All the sources mentioned in Sect. 18. li do not allow cross 
calibration down to the level of the sensitivity of future 
CMB experiments and do not enabl e us to detect sig- 
nals from the dark ages as predicted bv lBasu e t al. ( 200^ . 
Therefore here we want to discuss a method to cross cali- 
brate the frequency channel using the spectral distortions 
induced by the CMB dipole (see Sect.El). These spectral 
distortions can be predicted with very high accuracy on a 
level of ^ fiK. This is ~ 1000 times lower than the dipole 
signal and should therefore allow cross calibration to a 
very high precision. 

Every CMB experiment is measuring the fluctuations 
of the CMB intensity on the sky. Due to the nature of 
the physical processes producing these fluctuations in each 
point they are related to the fluctuations of the radiation 
temperature. In a map making procedure the measured 



nuctuations m the intensity are translated into the huctua- 
tions in the radiation temperature. Usually the aim of any 
map making procedure is to reduce statistic and system- 
atic error in order to obtain a clean signal from these CMB 
temperature fluctuations. As has been discussed in Sects. 
O the superposition of blackbodies with slightly different 
temperature in second order induces y-distortions in the 
inferred temperature differences. The amplitude of these 
distortions depends on the temperature difference and the 
chosen reference temperature. In Sects. Eland 13 it has been 
shown that due to the two basic observing strategies, abso- 
lute and differential measurements, the dipole can induce 
y-distortions with y-parameters up to y ~ 10"^. It was 
also shown in Sect. El that contributions from higher mul- 
tipoles are much smaller (y ^ 10^^^ — 10^^). 

The y-distortions induced by the dipole can contam- 
inate the maps produced in the map making procedure 
on a level higher than the sensitivity of the experiment. 
Therefore it is important to chose the map making pro- 
cedure such that spectral distortions are minimized. The 
discussion in Sects. 01 has shown, that for this purpose in 
absolute measurements of the CMB sky it is the best to 
choose the temperature of the reference blackbody close 
to the beam RJ temperature. This implies that for CMB 
experiment the best choice for the temperature of the in- 
ternal calibrator is the full sky mean temperature Tq. In 
order to minimize the spectral distortions arising due to 
the superposition of blackbodies in differential measure- 
ments the best choice for the reference temperature used 
to relate the intensity difference maps to the temperature 
difference maps is the RJ temperature of the combined 
temperature distribution of both beams. This optimal ref- 
erence temperature will be time dependent due to the var- 
ious scanning strategies and orientations of the spinning 
axis relative to the dipole during observations. 

The purpose of this Sect, is not to discuss details about 
map making procedures but to show, that there are ways 
to manipulate these CMB maps in order to make the spec- 
tral distortions arising due to the superposition of black- 
bodies become useful for cross calibration purposes. As the 
discussion in Sect. ^ has shown, for this issue it is better 
to move Ticf as far as possible away from Tq. The optimal 
choice of T-^ci depends on the sensitivity of the experiment 
and on the observing strategy: For absolute measurements 
Trof should be as close as possible to Tq in order to min- 
imize the induced y-distortions but on the other hand it 
should be chosen such that the induced y-distortion due 
to the CMB dipole is still measurable within the sensitiv- 
ity of the experiment. But in principle there is no strong 
constrain, since the induced y-distortions can again elimi- 
nated afterwards. For differential measurements it is possi- 
ble to chose the optimal reference temperature for calibra- 
tion issues independent of the best map making reference 
temperature, but here it is important to compare regions 
on the sky with maximal temperature difference, i.e. with 
maximal angular separation in the observed field. Below 
we now separately discuss CMB experiments using differ- 
ential measurements with full and partial sky coverage in 
more detail. 



iable 1. y-distortion: A J = lid-y m /iK lor y as given m 
the left column in some of the Planck spectral channels. 
Here yd = 2.6 • 10~'^, ?/opt,o = 3?/d and yopt = 12 yd- 
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Full sky CMB surveys 

For full sky missions like Planck the regions with max- 
imal temperature difference are located around the ex- 
trema of the CMB dipole. Using equation H75|l we see 
that they correspond to Afi = 2 and the resulting maxi- 
mal spectral distortion is characterized by yopt = 12 yd = 
3.1 • 10"^, ff we set Trof = Tr.i or Tj-of = i.e. to the 
maximum or minimum of the temperature on the CMB 
sky arising due to the dipole. 

One attractive procedure to compare the maximum 
and minimum of the CMB dipole is to take the CMB in- 
tensity maps of each spectral channel and to calculate the 
difference between each map and the map obtained by 
remapping the value of the intensity I(n) at position n 
to the value I{—n) at position — n, i.e rotating the initial 
map by 180 degrees around any axis, which is perpendic- 
ular to the CMB dipole axis and is crossing the origin or 
equivalently setting /i2 = — /ii in equation H67|l . 

Afterwards these artificial intensity maps for each 
spectral channel are converted into temperature maps 
using equation (0 and setting the reference tempera- 
ture at each point to r*(n) = T^^i{n) of the intrinsic 
map. Neglecting the contributions from intrinsic multi- 
poles with I > 2, the difference map in some frequency 
channel x^: = hv/kT-^^i will be given by 



Ar(x,,n) 



= -2/3dC+A^ + 2/3^C+/i' 

+ 12ydC+A^'-5y(i*), (78) 



where we used the definition jsj for gy(x). This shows, 
that the artificial maps will contain a frequency indepen- 
dent temperature dipole component with twice the ini- 
tial dipole amplitude and a monopole and quadrupole y- 
distortion resulting from the second term in jTHJ. The 
maxima of the spectral distortion will coincide with the 
extrema of the CMB dipole. 

Now, using the artificial map of the lowest frequency 
channel as a reference and taking the differences between 
this reference map and the artificial maps in higher spec- 
tral channels we can eliminate the frequency independent 
term corresponding to twice the dipole in equation (|78|) . 
This then opens a way to cross calibrate all the channels 
to very high precision, since the quadrupole component 
in these artificial maps will be larger in higher spectral 
channels, following the behavior of gy with x. 

The signal will include both statistical and system- 
atic errors for the average temperature of the sky and the 



dipoie amplitude p^. All uncertainties are mainly mnu- 
encing the frequency independent dipoie term but for the 
frequency dependent monopole and quadrupole they will 
become important only in next order: 

Ti./32=ri,o/3lo(l±AT)-(l±A^)2 

= Ti^o f3lo [1 ± 2 ± At] + 0{A^ ■ At) , (79) 

where Ti^o is the correct value of the temperature at a 
given point on the sky, /3d, o is the correct dipoie ampli- 
tude and At and A^ are their corresponding relative un- 
certainties. This implies that all the corrections to the fre- 
quency dependent terms are at least 1000 times smaller 
than the signal we are discussing here. On this level cor- 
rection due to higher order moments might become im- 
portant, but their contributions will be less than 1%. This 
means that y-distortion on a level of 10 /iK can be pre- 
dicted with ^ 100 nK precision. 

Let us note here, that applying the procedure as de- 
scribed above we are adding signals which have statis- 
tically independent noise. Therefore the statistical noise 
of the artificial maps should be a factor of 2/y/2 = y/2 
stronger than the initial maps. On the other hand we are 
increasing the amplitude of the y-distortion quadrupole 
by a factor of 12, resulting in strong gain. 

The maxima of the quadrupole component in the ar- 
tificial maps will be rather broad and correspond to thou- 
sands of square degrees on the sky, allowing to average 
the signal and thereby increasing the sensitivity. For cir- 
cular beam average formula (|78|l will be applicable. The 
statistical sensitivity of the Planck experiment will per- 
mit to find these spectral distortions in the artificial maps 
down to the level which is ne cessary to detect th e effects of 
reionization as discussed bv lBasu et alJ l)200 3l). Applying 
the above procedure to the 33 GHz and 93 GHz frequency 
channels of Wmap will lead to a maximal temperature dif- 
ference of AT = 3.6 fiK, whereas the maximal difference 
will be AT = 35.6 fiK for the 353 GHz channel of Planck 
(see TableOJ. Nevertheless, even in the case of Wmap the 
very high precision of the experiment might permit the 
detection of this quadrupole and thereby open the way to 
cross calibrate its different spectral channels. 

For Cobe/Firas the proposed method permits to 
check the precision of the internal calibration. The inter- 
nal calibrator was measured and tested with some finite 
precision and there is the possibility that the internal cal- 
ibration may be better than the precision at which it was 
tested. The proposed method has distinct spectral and an- 
gular properties, making it possible to improve the cross 
calibration of the spectral channels and possibly this will 
open a way to further improve the great results this ex- 
periment already gave us. 

CMB surveys with partial sky coverage 

Surveys with partial sky coverage may not simultaneously 
have access to regions around the maximum and the min- 
imum of the dipoie. Therefore here one should choose two 
regions in the accessible field of view such that the mean 
temperature difference between these regions is as large a 
possible. The dependence of the signal on the separation 
angle has been discussed in detail in Sect.[7| 



Beiore cross calibration all the bright sources detected 
such as bright clusters of galaxies and point sources in- 
side the chosen areas should be extracted. Typically the 
expected amount of SZ sources will be of the order of 
a few tens per square degrees with corresponding angu- 
lar extensions less than 1' (seelCarlstrom fc Holdeill2002t 
iRubino-Martin &: Sunvaevl 1200! ^). Afterwards the differ- 
ence of the signals in the two patches can be taken and 
compared to the signals in the other spectral channels. 

For example choosing one area centered on the min- 
imum or maximum of the CMB dipoie as a reference 
(T* = Ti-^i) and choosing the second area in the ring per- 
pendicular to the dipoie axis the maximal y-parameter is 
y = 3ydC+ = 7.7-10"'' C-i- (cf. equation (|77|l l. If we instead 
set — Tq we obtain the same maximal y-parameter but 
a stronger dependence on the beam radius or size of the 
regions we average (cf. equation (|75|l '). This small example 
shows that for any experiment with partial sky coverage 
a separate analysis of the optimal choice of the reference 
temperature and the regions on the sky has to be done. 



9. Other sources of spectral distortions 

Injection of energy into the CMB prior to recombination 
leads to spectral distortions of the background radiation. 
Given the Wmap best fit parameters, before redshift Zth ^ 
2 • 10^ all distortions are efficiently wiped out, whereas 
energy injection in the redshift range zth > z > z^, with 
z^ ^ 1.4 • 10^, leads to a /^-distortion and injection in 
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As was shown by ISilkl l|1968|) photon diffusion 
and thermal viscosity lead to the dissipation of 
small scale density perturbations before recombina- 
tion. This damping of acoustic waves will contribute 
to a fi- and y-dis tortion because it was lead i ng to 
(i) e n ergy release llSunvaev fc Zeldovichl Il970bt IPalvl . 
Il99lt iHu. Scott & Silkl Il994a|) and rii) mixing of 
photons from regions having differ ent temperatures 
(|Ze1dovich. Tllarionov Sun";^!???^. 

Now, the Wmap data implies that the initial spectrum 
of perturbations is close to a scale-invariant Harrison- 
Zeldovich spectrum w ith spectral index n = 1.03 ± 0.04 
l|Bennett et al.l l2003|) . Different estimates for the spectral 
distortions arising from the dissipation of acoustic waves 
in the early universe give a chemical potential of the order 
of Ur] \ s ~ 2 ■ 10"'' and y-distorti ons with j/dis ^ 10"'' l|Dalvl . 
ll99lllHu. Scott fc Silkl ^1994a^ . 

Another contribution to the y-distortion arise s from 
the epoch of re i onizati on l|Zeldovich fc Sunvaevl Il969t 
iHu. Scott fc Silkl . Il994b') and from the sum of the SZE of 
clusters of galaxies (Markcvitch et al, 1991). The Wmap 
results for the TE power spectrum point towards an early 
reionization of the universe with corresponding optical 
depth = 0.17 ± 0.04 (Kogut et al, 2003; SDcrgcl et3 
120031) . Here two main effects are important: (i) The pho- 
toionized gas typically has temperatures T of the order of 
~ 10"* K. Therefore the diffuse gas after reionization will 



produce a y-distortion m the CJVIB spectrum with 



kT 



dr - 3 • 10" 



T 



0.17 104 K 



(ii) The motion of the matter induces a y-distortion on the 
whole sky due to the secor id order Doppler effect, w ith cor- 
responding y-parameter ijHu. Scott fc SilkL Il994bj) . — 
i J (TT"-ec(i'^) dt, where {v'^) is the velocity dispersion 
over the whole sky in units of the speed of light. For a 
CDM model it is of the order of 
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2 • 10" 



10 



1/2 \ (v^) n^h^ 



- 1 



10-5 0.0224 



Since the reionization optical depth is very small, this is 
one order of magnitude smaller than j/re- If secondary ion- 
ization was patchy this will in addition lead to angular 
fluctuations ijSantos et al.U200^ . 

All the effects mentioned in this section lead to y- 
distortions with corresponding y-parameter in the range 
y ~ IQ-® — 10-*". As has been shown in Sects. and 
the y-distortions arising due to the CMB dipole are of the 
same order. For example the CMB dipole induces a full 
sky y-distortion with yd = 2.6 • 10"^. Therefore, to mea- 
sure any of the effects discussed above it is necessary to 
take spectral distortions associated with the CMB dipole 
into account. But since the angular distribution and the 
amplitude of these distortions can be accurately predicted 
they can be easily extracted. As has been shown in Sect. 
spectral distortions arising from higher multipoles typi- 
cally have y ~ 10~^^ — 10~^ and therefore do not play an 
important role in this context. But even for these distor- 
tions the locations and amplitudes can be accurately pre- 
dicted with measured CMB maps and therefore offer the 
possibility to eliminate these distortions from the maps. 

10. Conclusion 

We have discussed in detail the spectral distortions arising 
due to the superposition of blackbodies with different tem- 
peratures in the limit of small temperature fluctuations. 
The superposition leads to y-distortion with y-parameter 
ys = (S^) /2, where ((5^) denotes the second moment of 
the temperature distribution function, if the difference in 
temperatures in less than a few percent of the mean RJ 
temperature. We have shown that in this limit even com- 
paring to pure blackbodies leads to a y-distortion. 

The results of this derivation where then applied to 
measurements of the CMB temperature anisotropics with 
finite angular resolution and in particular to the CMB 
dipole and its associated spectral distortions, but in prin- 
ciple the method developed here can be applied when- 
ever one is dealing with the superposition of blackbod- 
ies with similar temperatures. We have shown, that tak- 
ing the difference of the CMB intensities in the direction 
of the maximum and the minimum of the CMB dipole 
leads to a y-distortion with y-parameter j/opt — 3.1 • 10-^. 
This value is 12 times higher than the y-type monopole, 
j/d = 2.6 • 10"^. Since the amplitude of this distortion can 
be calculated wi th the same precis i on as the CMB dipole, 
i.e. 0.3% today l|Fixsen fc Mather! . I2OO2I) . it opens a way 
to cross calibrate the different frequency channels of CMB 



experiment down to the level ol a lew tens ol nK ( lor more 
details see Sect.|Hl). 

We discussed another possibility to check the zero lev- 
els of different frequency channels by observing the differ- 
ence of the brightness in the direction of the dipole maxi- 
mum and in the direction perpendicular to the dipole axis: 
The dipole induced spectral distortion in this case is 4 
times weaker than for the difference of the maximum and 
minimum and corresponds to y ^ 7.7 • 10-^. Nevertheless, 
it is still 3 times stronger than the dipole induced whole 
sky y-distortion with y-parameter yd- 

The value of j/opt is only 5 times lower than the up- 
per limit on the whole sky y-parameter ob t ained by the 
Cobe/Firas experiment .Fixsen fc Mathe"!] (|200^ and is 
orders of magnitudes above the sensitivity of Planck and 
Cmbpol in each of their spectral channels. Therefore this 
signal might become useful for both cross calibration and 
even absolute calibration of different frequency channels 
of these experiments to very high precision, in order to 
permit detection of the small signals of r eionization as 
for example discussed bv lBasu et al l l|2003h . and to study 
frequency dependent foregrounds with much higher sensi- 
tivity. We should emphasize, that on this level we are only 
dealing with the distortions arising from the CMB dipole 
as a result of the comparison of Planck spectra with differ- 
ent temperatures close to the maximum and minimum of 
the CMB dipole, i.e. due to the superposition of blackbod- 
ies. The distortions are introduced due to the processing 
of the data and become most important in the high fre- 
quency channels. 

Let us stress that the distortions discussed in this work 
have both well known spectral and angular dependence 
and that they are connected with the much stronger signal 
of the CMB dipole. The amplitude of the dipole is known 
to a relative precision ~ 10-'^ and therefore the induced 
spectral distortions can be calculated to the same accu- 
racy. Maybe the great properties of the calibration source 
discussed above might even be detectable in the high fre- 
quency channels {x ^ 10) of the existing Cobe/Firas 
data and thereby help to further improve its calibration. 

The development of the CMB experimental technol- 
ogy is repeating with 20 years delay the progress that 
was made in the astronomical observations with X-ray 
grazing incidence mirrors and CCD detectors. Due to the 
big efforts of many spacecraft teams, the X-ray back- 
ground today is resolved to more than 85 % ijRosati et all 
.20021 . During the next decade, CMB observers will be 
able to pick up all rich clusters of galaxies and all bright 
y-distorted features connected with supernovae in the 
early universe, groups of galaxies and even rising down to 
patchy reionization in the CMB maps. Deep source counts 
may permit the separation of their contribution to the y- 
parameter and in principle might bring us close to the level 
of the y-distortions mentioned in Sect. El As has been ar- 
gued, on this stage the contributions to the y-parameter 
arising from the superposition of Planck spectra with dif- 
ferent temperatures corresponding to the observed CMB 
temperature ffuctuations will be the easiest to separate. 
Again, the main contribution will arise from the CMB 
dipole. On the full sky the corresponding y-parameter will 
be yd = 2.6 • 10"^. But even the distortions arising from 
the primordial temperature fluctuations with multipoles 
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t > 2 can lead to distortions ot the order or ?/ ~ 1(J~ 
in significant parts of the sky. In the future this type of 
spectral distortion can be taken into account making it 
possible to enter the era of high precision CMB spectral 
measurements. 



